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In this paper, we give random matrix theory approach to the quantum mechanics 
using the quantum Hamilton-Jacobi formalism. We show that the bound state prob¬ 
lems in quantum mechanics are analogous to solving Gaussian unitary ensemble of 
random matrix theory. This study helps in identify the potential appear in the joint 
probability distribution function in the random matrix theory as a super potential. 
This approach allows to extend the random matrix theory to the newly discovered 
exceptional polynomials. 



INTRODUCTION 


Q. 


Nine different formulations of non-relativistic quantum mechanics exist 111. These are 
the wavefunction, matrix, path integral, phase space, density matrix, second quantization, 
variational, pilot wave, and Quantum Hamilton-Jacobi (QHJ) formulations. It is interesting 
to note the close resemblance between random matrix theory and the QHJ formalism. The n 
stationary points of the random matrix ensemble action fl H = V(x)~lJ2k&log{\x k -xi\) 
where V(x) is the potential. In QHJ formalism, the problem is solved using the quantum 
momentum function (QMF) and it plays an important role . The quantum momentum 
function is defined as 3, 4] p — XT , -f 1 —b Q(x) here the moving poles are simple poles 
with residue —ih (we take here Ti — m — 1) and Q(x) is the singular part of the quantum 
momentum function. The moving poles are described as in the n first-order poles in the 
QMF of quantum mechanical origin. The positions of the poles are and same as the zeros 
of the wave function. As the zeros of the wave function change their positions with energy, 
so do the location of the corresponding poles in the QMF jdj]. The moving poles behave like 
random numbers as they are energy dependent (as the energy changes the position of these 
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moving poles will change). The above statement clearly captures the connection between 
the n stationary random points and the n moving poles of QMF. 

There are other models in the literature whose equations resemble the random matrix 
ensemble action. They are the Stieltjes electrostatic problem [5, 6] with n moving unit 
charges, interacting through a logarithmic potential, are placed between two fixed charges p 
and q at a and b respectively on a real line is given by V(xk) + Tk ~ Xl = 0 where V(xk) 
is the position of the two fixed charges. The motion of point vortices in hydrodynamics 
is governed by the Helmholtzs equations|7] = A- i + W(zk) here the R, is 
circulations or strengths of point vortices at Zj respectively and W(zk) is the back ground 
flow. The equation for low-energy spectrum of the Heisenberg ferromagnetic chain [8] ^ = 
nd+Y Yjm^l ui-Ur, -. The analogy between the Stieltjes electrostatic problem and the quantum 
momentum function is established by the author him self 9] and the connection between the 
deisenberg ferromagnetic chain and the Stieltjes electrostatic problem was established in ref 

B|. 


In this paper, we show that the solving quantum mechanics using QHJ formalism is 
analogous to solving the quantum mechanics problem using the random matrix theory. The 
paper is organised as follows, in the following section of the introduction we briefly describe 
the RMT (part A) followed by QHJ (part B) and Supersymmetric quantum mechanics 
(part C). In the section II, we give the equivalence between the QHJ and RMT followed by 
conclusion in section III. 


A. Random Matrix Theory 


In 1951, Wigner proposed the use of random matrix theory (RMT) to describe the prop¬ 
erties of excited states of atomic nuclei 0 . This was the first time RMT was used to model 
physical reality. In random matrix theory, the dynamics of the ensemble of a random infinite 
dimensional Hermitian matrix is described by the probability distribution function 


V(Xi ,..., Ajv) dA — c n e 1311 dA 


( 1 ) 


where dA = dXi... dX tv, c n is constant of proportionality and the index (3 = 1,2,4 char¬ 
acterizes the real parameters of the symmetry class of orthogonal, unitary and symplectic 
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respectively for the random matrix 

N N 

# =-5>(Ai)-5> (2) 

i= 1 i<l 

As the random matrix is invariant under the symmetry group these ensembles are called the 
invariant ensembles. In the random matrix H, V(\) is the potential, A,; are the eigenvalues 
with i being a free index running from 1, 2, • ■ • N and \X t — Xj\ is the Vandermonde determi¬ 
nant. For details please refer to [3]. The probability distribution function can be rewritten 
as 

V(Xi, ...,X N )dA = c n e~^=^ v ^ I] l A * - A / dA ( 3 ) 

i<j 

By manipulating the Vandermonde determinant, that is, by adding and deleting columns or 
rows, the equation (J3J) is written as 

N 

■p(Ai,..., A at) dA = c n \\_ w^(Xi) n | A, - Aj | dA (4) 

i— 1 %<j 

where wp(Xi) is the weight function of classical orthogonal polynomials. The classical orthog¬ 
onal polynomials are classified into three different categories depending upon the range of 
the polynomials. The polynomials in the intervals (— oo; oo) with weight function w = e~ x2 
are the Hermite polynomials. In the intervals [0; oo) with weight function w = X b e~ x are the 
Laguerre polynomials. In the intervals [—1; 1] with weight function w = (1 + A) a (l — X) b are 
the Jacobi polynomials. 

The invariant ensembles in random matrix theory are classified into three the Gaussian 
ensembles, the Wishart ensembles and the two Wishart ensembles. In the case of Gaussian, 
the matrices are known as the Wigner matrices. The Gaussian ensembles are further classified 
the Gaussian orthogonal ensemble (GOE), the Gaussian unitary ensemble (GUE), and the 
Gaussian symplectic ensemble (GSE). In stationary case, for the Gaussian ensembles, Wigner 
law states that the empirical distribution of eigenvalues Ai,..., Aat as N tends to infinity 
converges to a semicircle in the interval [a,b]. In non-stationary case the equation (J3]) is the 
stationary solution of a Fokker-Planck equation Q. 

In this paper we only deal with the Gaussian ensembles and show that the bound state 
problems in quantum mechanics are analogous to solving Gaussian unitary ensemble of 
random matrix theory. This study helps in identifying the potential that appear in the joint 
probability distribution function in the random matrix theory as a super potential. This 
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approach allows to extend the random matrix theory to the newly discovered exceptional 
polynomials. 


B. Quantum Hamilton Jacobi 


In this section, a brief review of Quantum Hamilton Jacobi formalism is presented below. 


For details see the references 

ti 2 


13]. The Schrodinger equation is given by, 


- ^V 2 ^(x, y-> z ) + V ( x > Vi z )^i x , V, z ) = y, z ). 


(5) 


One defines a function S analogous to the classical characteristic function by the relation 


ip(x,y,z) =exp (S.\ 


( 6 ) 


which, when substituted in (f5j), gives 

(VS) 2 -ihV.(VS) = 2 m(E-V(x,y,z)). 
the quantum momentum function p is dehned in terms of the function S as 

p= VS. 

Substituting (jHJ) in ([7j) gives the QHJ equation for p as 

ip) 2 ~ ihV.p = 2 m(E — V (x, y, z)) 

and from ([5]) and (JHJ) , one can see that p is the the logarithmic derivative of -0(x, y, z) i. e. 

p = —ihVlniJj{x,y, z) (10) 

The above discussion of the QHJ formalism is done in three dimensions the same equation 
in one dimension takes the following form 


(7) 


( 8 ) 


(9) 


p 2 — ift—j— = 2 m(E — V(x)), 
ax 


( 11 ) 


which is also known as the Riccati equation. In one dimension the eq (flUjl take the form 


p = —i% — Inil’(x). 
ax 


( 12 ) 

ax 

It is shown by Leacock and Padgett jll, 12] that the action angle variable gives rise to exact 
quantization condition 

J{E) = — i pdx = nh. 

Zn Jc 


( 13 ) 






5 


C. SUSY Quantum mechanics 


In supersymmetry, the superpotential W(x) is defined in terms of the intertwining oper¬ 
ators A and A 1 as 

A=4- + W(x), A^ = ~ + W(x), (14) 

ax ax 

This allows one to define a pair of factorized Hamiltonians H ± as 

d 2 


H + = A'A = 
H~ = Hit = 


dx 2 

d*_ 

dx 2 


+ V + (x) -E, 
+ V~{x)-E, 


(15) 

(16) 


where E is the factorization energy. The partner potentials V ± (x) are related to W(x) by 

V ± (x) = W 2 (x) =F W{x) + E, (17) 


where prime denotes differentiation with respect to x,refer the reader to 15] for more details. 


II. EQUIVALENCE BETWEEN RMT AND QHJ 


To establish the relation between the RMT and QHJ we start with the quantum momen¬ 
tum function given in equation (11111 . In equation ([Till p is the QMF and is given by equation 
(lT2j) . It is evident form the equation (1T2]) the QMF is related to wave function. Further 
more the relation between the QMF and the wave function in the equation (fT2j) captures the 


randomness arising from the moving poles in QMF. 


in terms of n moving poles and fixed poles 


fly, 
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hen the quantum momentum function 


14] is given by 


n —if) 
p(x) = £ 

mX-Xfc 


Q(x). 


By introducing the polynomial 


/(x) = (x - xQ(x - x 2 ) • • • (x - x n ), 

then the quantum momentum function ( h — 1) in terms of the polynomial reads as 


(18) 


(19) 


( 20 ) 


and substituting in (fill) then one gets 


fix) + 2 iQ(x)f\x) + [Q 2 {x) - iQ'(x) - E + V(x)]f(x) = 0. 


( 21 ) 
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The search for the polynomial solutions to the equation (12Tj) leads to quantization. This 
is equivalent to demanding [ Q 2 (x ) — iQ'{x ) — E + V{x)\ to be constant. This will only be 
possible if Q(x ) = zW(x). In QHJ it turns out that the Q(x) is the super potential. It 
should be pointed out here that the energy in the equation (12T1) id different from the energy 
in equation (ITTjh Then the equation (T2T1) reduces to 

- f"(x) + 2 W(x)f'(x) = —E n f(x). (22) 


It is evident from the above equation that the left hand side ( Ihs) balances the right hand 
side ( rhs ). By dividing equation (1221 ) with 2/'(x fc ) and defining the Ihs of equation (12211 as 
T(x) one gets 


Tix) - J M +w{x) ' 


(23) 


Now, we will show that QHJ is analogous to the Steiljes electrostatic model and then 
show that QHJ is equivalent to RMT. By following the argument in the Metha’s book 

1 f(x) 1 (x - Xj)f\x) - f(x) 


E 


x - x k f{x) X ~ Xj (X - Xj)f(x) 

by taking the following limit x —>■ Xj and using l’Hospital rule one gets 


(24) 


E 


1 <j<n,j^k X 3 Xk 


= lim 

X^>X 


f( x ) 


= lim 

x—>x. 


j lf(x) X — Xj 
(x - Xj)f'(x) - f(x) 


(x - Xj)f(x) 


nxj) 


(25) 


Substituting 


m 


one gets 


T{x) = E 


W(x), 


_ _ ... (26) 
1 <j<n,j¥=k Xk X 3 

which is minimum of the Lagrangian considered for the the Steiljes electrostatic model 0 . 0 . 

By re-substituting Q(x) = iW(x ) in equation (1261) the quantum momentum function 
reduces to 


p( x j) = E 


—i 


+ Q(xj). k = 1, 2 • • ■ n. 


(27) 


1 <j<n,j^k Xk X 3 

In order to compare the RMT and QHJ, we define the probability distribution function 
in RMT as the wave function 


'/'(V) = P(Ai, ...,\ N )= C n exp BH 


(28) 
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where 

H=V( AJ-i £ ("(l^-Ajl) (29) 

Z 1 <j<n,j^k 

Substituting (1281) in (!T2l) the quantum momentum function is given by 

p( A fc )= E -'S'(Afc) k = l,2---n (30) 

where S(xk) = ■ By comparing equation (]30l) and (lT8]i one gets (3 = — ih and W(x) = 

jjS(x) = Q(x). It is therefore clear that the potential used in the random matrix theory is 
the super potential W(x). Substituting the equation (TT9|) in equation (130]) and the substuting 
this in equation (fill) one gets the following differential equation as 

f"(x) + 2 iQ(x)f(x) + [Q 2 {x) - iQ'(x) - E + V(x)]f(x) = 0. (31) 


By exploting the supersymmetric quantum mechanics one gets the solutions for the above 


differential equation as the classical ort 


potential. For the Harmonic oscillator 
the differential equation (13TT) becomes 


;onal polynomials for the given range fixed by the 
15] the super potential is given by Q(x) = \ojx then 


/"(C) - 2C fix) + 2n/(C) = 0 1 < k < n. 


(32) 


where ( = ax and a 2 = The differential equation (IT?!) is the Hcrmite differential equation 
and the solutions are the Hermite polynomials given by ip(x) = exp[— \x 2 ]H n (x). One can 
show by supersymmetric arguments the potential is given by V(x) = \ux 2 . This model in 
RMT is known as the Dyson gas model. 

One can extend the RMT models to the recently discovered exceptional polynomials, as 
quantum Hamiltonian Jacobi formalism for the deformed oscillator is studied in ref [bij and 
the momentum function is given by 


P= E 


— x — 


= ] X ~ x k 


X — X\ 


+ 


g + l 


X 


(33) 


and thus the solutions are exceptional Laguerre polynomials we obtain the solution as ex¬ 
ceptional polynomials 


= 


X 9+l e 2 + 


-9+l-i 


(x 2 ) 


L? 


(x 2 ,g)- 


(34) 












It is clear from quantum mechanic the wave function is given by ip = w^(x)Pi(x) where 
w(x) is the weight function for the corresponding orthogonal polynomial P/. The probability 
of finding the particle in an interval is given by 

p — f ip(x)*ip(x)dx = I w(x)Pi{x)P m {x)dx. (35) 

Ja Ja 

In random matrix theory this turns out to be the Joint probability distribution function 
Gaussian unitary ensemble in terms of eigenvalues 

P(Ai,..., Ajv) dA — C n w\A(x)\ 2 dA. (36) 


where |A(x)| = T\, <3 |Aj — Aj| is the Vandermonde determinant. This the Joint probability 
distribution function for exceptional Laguerre polynomials in terms of eigenvalues is given 
by 


^(Ai,..., \n) dA 


r 2 (9+0 P -P 

Cn 2{g+l) _ 3 n l A i - A J '| 2 clA 

{X ) %< j 


(37) 


where L”(x 2 ,^) are exceptional Laguerre polynomials. 


III. CONCLUSION 

In this paper, we given random matrix theory approach to the quantum mechanics using 
the quantum Hamilton-Jacobi formalism. We shown that the bound state problems in 
quantum mechanics are analogous to solving Gaussian unitary ensemble of random matrix 
theory. We have also shown that the potential appear in the joint probability distribution 
function in the random matrix theory as a super potential. Using this approach we have 
extend the random matrix theory to the newly discovered exceptional polynomials. 
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